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1.  Introduction 


Depth-sensing  indentation,  also  referred  to  as  nanoindentation,  is  a  technique  in  which  a  hard 
indenter  tool  (usually  made  of  diamond  or  crystalline  AI2O3)  is  pressed  into  a  sample  surface 
while  load  and  displacement  are  measured  continuously.  It  is  commonly  employed  to  find 
elastic  modulus  and  indentation  hardness  properties  of  various  materials  on  the  microscale 
(1-3).  Because  indentation  crudely  resembles  a  uniaxial  compression  test,  it  may  seem  that 
stress-strain  relationships  may  be  readily  discerned  from  load-displacement  data  for  any  tool 
geometry.  However,  it  has  been  shown  that  such  a  relationship  is  unique  to  nonself-similar 
indenter  geometries  (i.e.,  a  sphere)  but  not  to  pyramidal  or  conical  indenters  (4). 

Pioneering  work  to  determine  expressions  for  indentation  stress  and  strain  was  done  by  Tabor 
and  focused  primarily  on  the  response  of  metals  to  contact  loading  (5).  Since  that  time,  a 
number  of  studies  have  attempted  to  describe  indentation  stress-strain  curves  with  some 
promising  success  (6-12).  These  works  have  concentrated  on  a  limited  modulus  range  or 
specific  class  of  material.  However,  no  general  study  exists  that  proposes  an  easily 
implemented,  general  framework  to  estimate  the  elastic  modulus  and  yield  stress  of  a  material 
for  spherical  indentation,  exclusively  using  the  loading-curve  data.  Such  a  procedure  is 
developed  in  this  work,  and  a  number  of  materials  important  for  military  applications  are  used  to 
exemplify  the  approach.  This  work  is  especially  useful  in  measuring  mechanical  stress-strain 
curves  on  materials  that  cannot  be  measured  through  bulk  measurement  or  materials  that  are 
heterogeneous  at  the  microscale.  Examples  include  single  material  grains,  thin  films,  and 
microelectromechanical  system  components. 


2.  Theory 


To  derive  the  relationship  for  indentation  stress  and  strain,  the  existing  contact  mechanics 
framework  is  considered.  The  first  analysis  of  a  spherical  body  orthogonally  indenting  an  elastic 
halfspace  was  done  by  Hertz  (13).  Figure  1  shows  a  schematic  of  many  variables  used  in 
equations  1-5.  The  elastic  displacement  of  a  sphere  into  a  material  may  be  expressed  by 


h  = 


3P 

V4  E  j 


\Rj 


(1) 
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Figure  1.  Schematic  of  indenter  contact  with  a  sample  surface.  (A  number  of  variables  from  equations 
1-5  are  depicted.) 


* 

where  P  is  the  load  applied  to  the  indenter,  R  is  the  sphere  radius,  and  E  is  the  reduced  modulus 
of  the  material,  defined  as 

E  E 

l 

where  vand  E  are  the  Poisson’s  ratio  and  elastic  modulus  of  the  indented  material,  respectively, 
and  the  subscript  i  refers  to  indenter  properties  (for  diamond,  Et  is  commonly  taken  as  1141  GPa 
and  v;.  is  commonly  taken  as  0.07). 


(2) 


For  purely  elastic  contact,  it  has  been  recognized  (14)  that  if  the  spherical  shape  is  approximated 
by  a  paraboloid  of  revolution,  the  contact  depth  can  be  related  to  total  elastic  penetration  depth 
by 


h,  = 


(3) 


As  soon  as  elastic-plastic  behavior  develops,  there  exists  a  plastic  depth  hp,  in  addition  to  the 
elastic  penetration  depth.  Therefore,  the  total  penetration  depth  becomes 

h,  =  K  +  hp  .  (4) 
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Because  the  plastic  depth  offers  no  component  in  the  elastic  deflection  of  the  material,  the  total 
contact  depth  will  be 


h 


c 


(5) 


The  contact  depth  may  be  determined  by  the  commonly  used  relation 


h 


c 


(6) 


where  fis  a  geometric  constant  that  is  0.75  for  a  paraboloid  of  revolution  (14)  and  S  is  the 
contact  stiffness  measured  continuously  throughout  the  duration  of  loading  by  the  continuous- 
stiffness  technique  (75).  The  continuous  stiffness  is  determined  by  superimposing  a  small 
oscillating  excitation  force  amplitude  F0  at  a  frequency  coo  n  top  of  the  quasistatic  applied  load 
and  using  the  relation 


-i 


1 


1 


^cos  </)-(Ks-Mco2)  Kf 
z„ 


(V) 


where  zQ  is  the  displacement  amplitude  resulting  from  F0,  (j)  is  the  phase  angle  shift  of  (O,  Ks  is 
the  support  spring  stiffness,  Mis  the  mass  of  the  indenter  column,  and  if/- is  the  frame  stiffness. 
Values  for  Ks,  M,  and  Kf  are  determined  by  the  response  of  a  free-hanging  indenter  tip.  The 
other  variables  in  equation  7  are  directly  measured  continuously  during  the  test,  so  hc  in 
equation  6  becomes  an  experimentally  determinable  value,  since  ht  and  P  are  also  directly 
measured. 


By  combining  equations  4  and  5,  and  substituting  in  for  he  via  equation  1,  the  relation 


K-h 


c 


3  P 

4  E* 


,  2/3 


f  i  Y/3 


(8) 


is  obtained,  with  any  displacement  due  to  plasticity  inherently  accounted  for.  Note  that  this 
equation  reduces  to  equation  1  when  displacement  is  purely  elastic.  One  fact  often  ignored  in 
indentation  analysis  is  that  the  radius  of  spherical  and  rounded  conical  indenter  tips  is  often  not 
constant  nor  is  it  the  quoted  nominal  value  (16).  Therefore,  it  is  imperative  to  allow  the  value  of 
R  to  vary  as  a  function  of  depth  as  it  appears  in  equation  8  to  more  adequately  describe  real 
indenters.  It  can  be  shown  geometrically  that  R  may  be  related  to  the  contact  depth  and  projected 
contact  radius  a  by 


R  = 


a2  +hc2 

2 K 


(9) 
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For  any  indenter,  a  will  be  a  function  of  hc,  related  to  the  calibrated  area  coefficients  Cn. 


Combining  equations  2,  8,  and  9  and  solving  for  the  elastic  modulus  of  the  material, 


4 

2(h,-hc) 

3/2 

1  2  .  (11) 

1-1/,.  V  7 

3  P 

(  2  hr  ] 

1/3 

E, 

2  .  r  2 
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By  definition,  the  elastic  modulus  is  the  ratio  of  the  true  stress  over  the  true  strain.  As  long  as 
true  stress  can  be  determined,  then  an  expression  for  true  strain  will  follow.  The  definition  of 
true  indentation  stress  is 


a  = 


y/P 

— T  ’ 

7ta 


(12) 


or  load  over  projected  contact  area,  where  (//is  related  to  a  parameter  called  the  constraint  factor 
(5).  This  factor  has  been  found  to  be  ~0.9  on  average  when  a  material  is  in  the  purely  elastic 
regime,  decreasing  to  0.61  when  in  the  elastic-plastic  regime,  and,  finally,  tending  to  0.35  when 
full  plasticity  has  been  developed  ( 1 7).  In  the  present  work,  it  has  been  found  experimentally 
that  a  suitable  expression  to  describe  the  behavior  of  ^through  the  three  loading  regimes  is 


y/  =  0.9 


\K  j 


(13) 


The  ratio  of  h/hc  will  be  equal  to  2  if  the  contact  is  purely  elastic  and  equal  to  1  if  the  contact  is 
purely  plastic.  The  value  of  t//  will  tend  to  zero  if  a  purely  plastic  condition  arises,  unlike  tending 
to  0.35,  as  in  Francis  (77).  This  was  done  for  simplicity  and  because  it  more  accurately 
described  the  results.  Because  the  indentation  stress  decreases  significantly  when  plastic 
behavior  is  the  dominating  mechanism,  equation  13  greatly  diminishes  in  its  accuracy  when 
testing  extremely  plastic  materials  with  a  relatively  small  sphere,  e.g.,  aluminum  with  a  3 -pm 
nominal  radius  indenter.  When  equation  13  is  substituted  into  equation  12,  the  expression  for 
true  indentation  stress,  which  is  used  throughout  this  work,  is 


(7  = 


0.9  P(h,  ') 

- 2~ 

m  VK  j 


(14) 
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Recognizing  that  the  factor  ^should  be  linked  with  P  and  not  projected  contact  area,  the 
expression  for  true  strain  used  throughout  this  work  becomes 


s  = 


f  0.9 P  ihL_  ^ 
m2(l-v2)\hc 


2.7  P 


^-1 

\K  . 


2{h,-he) 

/  \  1/3 

2  h  A 


va“  +  V  y 


3/2 


1-v. 


(15) 


All  variables  in  this  expression  depend  only  on  measurements  of  load,  total  displacement  into 
surface,  and  harmonic  contact  stiffness. 


3.  Experimental 


Compression  and  indentation  tests  were  done  on  cylindrical  faces  of  five  different  materials  with 
the  following  dimensions:  ~10-mm  height  and  5-mm  diameter,  yielding  a  length-to-diameter 
ratio  of  2: 1 .  The  materials  used  for  study  were  two  polymers  (polycarbonate  [PC]  [Lexan  9034 
from  GE  Plastics]  and  polymethylmethacrylate  [PMMA]  [Plexiglas  G  from  Atofina  Chemicals]), 
two  metal  alloys  (rolled-homogeneous-armor  [RHA]  steel  [75]  and  a  titanium-aluminum6%- 
vanadium4%  [Ti-6-4]  alloy),  and  a  ceramic  composite  made  from  tungsten  carbide  and  12%-by- 
volume  cobalt  (WC-Co).  These  materials  were  chosen  because  of  their  range  in  mechanical 
properties  and  widespread  use  in  U.S.  Army  applications. 

Uniaxial  compression  tests  for  all  materials  were  perfonned  with  true  strain  controlled 
compression  on  an  Instron  1331  (MTS  Systems).  This  mode  was  achieved  by  using  a 
programmable  function  generator  that  sent  the  command  signal  feed  to  the  auxiliary  input  of  the 
Instron  machine.  Liberal  amounts  of  a  lubricant  were  applied  to  each  specimen  to  reduce  friction 
between  the  specimen  ends  and  platens.  Also,  small  grooves  were  made  on  the  faces  of  the 
Ti-6-4  and  RHA  steel  specimens  to  retain  the  lubricant  to  the  greatest  compression  strains 
possible.  No  volume  change  was  assumed  throughout  the  end  of  the  tests  under  pressure.  Tests 
were  perfonned  at  a  strain  rate  of  0.001  1/s.  Compression  tests  frequently  show  evidence  of 
nonperfect  contact  at  the  beginning  of  the  curve  data  for  most  specimens,  so  only  compression- 
yield-stress  values  were  compared  with  indentation  data.  Elastic-modulus  values  for  each 
material  were  detennined  by  ultrasonic  techniques,  as  was  the  case  for  the  Ti-6-4,  RHA,  and 
WC-12Co,  or  by  dynamic  mechanical  analysis  (DMA)  at  0.1  Hz,  as  was  the  case  for  the  PC  and 
PMMA  samples. 
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Indentation  tests  were  run  using  a  Nano  Indenter*  XP  (MTS  Systems),  which  has  a  load  sensing 
resolution  in  the  fiN  range  and  depth  sensing  resolution  in  the  sub-nm  range.  Before  indentation, 
all  samples  were  mounted  in  epoxy  and  carefully  polished  to  a  mirror  finish  as  seen  under  a  50x 
microscope  objective,  using  0.5-p.m  diamond  particle  suspensions  in  the  final  polishing  step.  A 
Dimension  3100  atomic  force  microscope  (Digital  Instruments)  was  used  in  tapping  mode  to 
determine  the  average  root-mean-square  roughness  value  for  each  sample  over  a  20-  x  20-pm 
scan  area.  To  eliminate  compliance  effects  due  to  the  epoxy  mount,  samples  were  ground  from 
the  backside  until  the  sample  material  was  exposed  and,  finally,  were  adhered  to  a  stiff 
aluminum  puck.  All  tests  were  run  at  an  indentation  strain  rate  (loading  rate  divided  by  load)  of 
0.05  1/s,  and  the  maximum  allowed  thermal  or  electronic  drift  rate  was  0.05  nm/s.  Different 
strain  rates  were  attempted  but  either  not  enough  data  points  were  collected  at  shallow  depths  or 
drift  effects  would  begin  to  significantly  skew  the  results.  Three  different  indenter  tips  were 
used:  two  diamond  sphero-conical  geometries  with  20-  and  50-p.m  nominal  radii  and  one  ruby- 
spherical  geometry  with  a  500-pm  radius.  For  each  geometry  and  sample,  a  4  x  4  array  of 
indents  spaced  at  least  100  |im  apart  was  generated. 


4.  Results  and  Discussion 


In  table  1,  the  results  from  bulk  mechanical  testing  (compression,  ultrasonic,  and  DMA  testing) 
are  summarized  in  terms  of  compressive  yield  stress  and  elastic  modulus  values.  In  tables  2M, 
estimated  compressive  yield  stress  and  elastic  modulus  values  are  shown  for  instrumented 
indentation  measurements  using  20-,  50-,  and  500-pm-radius  tips.  Compressive  yield  stress 
values  were  determined  from  compression  and  indentation  testing  such  that  they  represent 
(approximately,  in  the  case  of  indentation)  offset  yield  stress  values  at  a  strain  offset  of  2%.  The 
compression  values  shown  in  table  1  were  measured  and  calculated  using  standard  methods  and, 
thus,  are  estimates  of  the  actual  yield  stress  values  to  which  the  indentation  results  are  compared. 

For  the  two  softer  materials,  PMMA  and  PC,  elastic  modulus  values  determined  from 
indentation  for  each  of  the  three  indentation  tips  (tables  2-4)  were  in  good  agreement  with  the 
elastic  modulus  measured  using  DMA  (table  1).  While  this  agreement  is  expected  of  the 
continuous  stiffness  measurement  (CSM)  method,  it  shows  that  the  indentation  stress-strain 
curves,  constructed  from  the  indentation  force-displacement  data  using  equations  14  and  15, 
have  the  correct  slope.  For  PMMA  and  PC,  modulus  values  determined  using  the  20-pm-radius 
tip  were  closest,  in  general,  to  the  DMA  measurements.  However,  for  the  harder  materials 
(Ti-6-4,  RHA  steel,  and  WC-Co),  the  elastic  modulus  values  from  the  CSM  method  and 
indentation  stress-strain  slope  method  were  closest  to  the  ultrasonic  method  (table  1)  when  using 


Nano  Indenter  is  a  registered  trademark  of  MTS  Systems  Corporation. 
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Table  1.  Compressive  yield  stress  values  for  Ti-6-4,  RHA  steel,  WC-Co,  PMMA,  and  PC  measured  using 

compression  testing  and  elastic  modulus  values  measured  using  ultrasonic  testing  (Ti-6-4,  RHA  steel, 
and  WC-Co)  and  DMA  (PMMA  and  PC). 


Material 

Ti-6-4 

RHA  Steel 

WC-Co 

PMMA 

PC 

Compression  Yield  (GPa) 

1.38 

1.22 

4.51 

0.113 

0.0682 

Bulk  Measured  Modulus  (GPa) 

114 

205 

530 

E'  =  4.24  E"  =  0.323 

E'  =  2.21  E"  =0.029 

Table  2.  Elastic  modulus  and  compressive  yield  stress  values  estimated  from  indentation  testing  of  Ti-6-4,  RHA 
steel,  WC-Co,  PMMA,  and  PC  using  a  20-gm-radius  indentation  tip. 


Material 

Ti-6-4 

RHA  Steel 

WC-Co 

PMMA 

PC 

Compression  Yield  (GPa) 

1. 1-1.2 

0.9-1. 2 

5.5-6 

0.065-0.075 

0.06-0.08 

Stress-Strain  Curve  Modulus  (GPa) 

65-70 

160-180 

330-350 

2.5-3 

1.9-2 

CSM  Average  Modulus  (GPa) 

97  ±2 

209  +  2 

415  +  15 

4  +  0.4 

2  +  0.5 

Table  3.  Elastic  modulus  and  compressive  yield  stress  values  estimated  from  indentation  testing  of  Ti-6-4,  RHA 
steel,  WC-Co,  PMMA,  and  PC  using  a  50-pm-radius  indentation  tip. 


Material 

Ti-6-4 

RHA  Steel 

WC-Co 

PMMA 

PC 

Compression  Yield  (GPa) 

0.9-1. 3 

1-1.1 

4.5-5 

0.03-0.04 

0.05-0.06 

Stress-Strain  Curve  Modulus  (GPa) 

85-95 

180-185 

420-440 

2.7-2. 8 

1. 8-2.5 

CSM  Average  Modulus  (GPa) 

110  +  3 

232  +  5 

520  ±  19 

3.4  ±  1.3 

1. 7  +  0.1 

Table  4.  Elastic  modulus  and  compressive  yield  stress  values  estimated  from  indentation  testing  of  Ti-6-4,  RHA 
steel,  WC-Co,  PMMA,  and  PC  using  a  500-gm-radius  indentation  tip. 


Material 

Ti-6-4 

RHA  Steel 

WC-Co 

PMMA 

PC 

Compression  Yield  (GPa) 

Unmeasurable 

Unmeasurable 

Unmeasurable 

Unmeasurable 

Unmeasurable 

Stress-Strain  Curve  Modulus  (GPa) 

40-50 

— 

350-400 

2.5-2. 6 

1.3-2 

CSM  Average  Modulus  (GPa) 

55  ±3 

— 

381+60 

5.5  ±0.3 

2.2 +  0.1 

the  50-pm-radius  tip  (table  3).  With  the  exception  of  the  measurements  on  PC,  the  greatest 
deviation  in  modulus  values  between  indentation  and  bulk  mechanical  measurements  was 
observed  when  using  the  CSM  technique  with  the  500-|nn-radius  tip. 

Compressive  yield  stress  values  for  indentation  and  bulk  measurement  were  in  general 
agreement  with  one  another.  For  almost  all  materials,  there  was  less  than  33%  difference 
between  the  two  measurement  techniques  (except  for  PMMA,  where  the  difference  was  as  much 
as  73%).  The  20-  and  50-pm-radius  indenters  estimated  the  compressive  yield  stress  effectively. 
For  PMMA  and  PC,  the  20-pm-radius  tip  gave  values  closest  to  the  compression  tests.  For 
WC-Co,  the  50-pm-radius  indenter  gave  the  best  values.  In  all  cases,  the  500-pm-radius 
indenter  was  not  successful  at  capturing  a  compressive  yield  stress  value  because  of  the  low 
overall  stress  (as  described  by  equation  13)  at  maximum  load.  For  this  reason,  it  is  advised  to 
use  spheres  with  radius  values  in  the  range  of  10-100  pm  in  order  to  reach  a  suitable  stress 
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state  underneath  the  tool  that  would  make  the  analysis  applicable.  A  3-um-radius  tool  was  used 
to  get  compressive  yield  stress  values  via  this  method  (data  not  shown),  but  results  were 
inconclusive  and  very  scattered.  Therefore,  it  was  found  that  if  the  stress  generated  underneath 
the  indenter  is  too  little  or  too  great  at  shallow  depths,  this  method  does  not  provide  good  results. 
This  conclusion  pertains  to  soft  and  hard  materials. 

Surface  roughness  played  a  large  role  in  the  accuracy  of  results.  Although  all  samples  were 
meticulously  polished  to  achieve  the  smoothest  surface  possible,  roughness  values  were  different 
between  the  samples.  Average  root-mean-square  roughness  values  measured  were  2.4  nm  for 
Ti-6-4,  3.9  nm  for  PC,  7.6  nm  for  RHA  steel,  12.8  nm  for  WC-Co,  and  16.6  nm  for  PMMA. 
PMMA,  the  sample  with  highest  surface  roughness,  was  the  most  difficult  for  indentation  to 
describe  accurately.  For  all  other  samples,  measurements  with  the  the  20-  or  50-um-radius  tip 
overlapped  compression-test  results.  When  there  was  high  surface  roughness  for  a  given  sample, 
it  was  expected  that  the  modulus  and  compressive  yield  stress  would  be  lower  than  expected. 
This  is  because  surface  asperities  create  a  nonperfect  contact  between  the  tip  and  sample, 
increasing  the  contact  compliance  and  lowering  modulus  and  compressive  yield  stress.  In  fact, 
except  for  the  compressive  yield  stress  of  WC-Co  with  the  20-pm  sphere,  all  loading-curve 
modulus  and  yield-stress  values  were  lower  for  indentation  compared  to  compression.  This 
effect  was  likely  attributable  to  the  asperities. 

Data  showing  stress-strain  behavior  as  measured  by  compression  testing  with  the  20-um-radius 
indenter  is  shown  in  figures  2a-e.  Again,  because  of  imperfect  initial  loading  contact,  the 
compression  data  was  not  used  to  calculate  modulus,  only  compressive  yield  stress.  However, 
indentation  was  able  to  capture  both  measurements.  In  all  cases,  it  was  seen  that  the  highest 
indentation  stress  reached  was  similar  to  that  of  the  compression  data.  Similar  curves  were 
generated  for  the  50-um-radius  indenter.  With  the  500-um-radius  indenter,  curves  were  also 
similar  except  for  no  clear  maximum  yield  stress  in  all  cases.  For  all  indenters,  results  were 
sensitive  to  where  the  initial  surface  contact  (h  =  0  mn)  was  taken  to  be.  This  point  was  determined 
automatically  via  a  defined  hannonic  contact  stiffness  value  set  by  the  Nano  Indenter  XP. 
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True  Strain  (mm/mm) 

(a) 


True  Strain  (mm/mm) 

(b) 


Figure  2.  Curves  showing  true  stress  and  true  strain  for  compression  data  (solid 
lines)  and  20-pm-radius  indenter  data  (dotted  lines).  Curves  are  shown 
for  (a)  Ti-6-4,  (b)  RHA  steel,  (c)  WC-Co,  (d)  PMMA,  and  (e)  PC. 
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Figure  2.  Curves  showing  true  stress  and  true  strain  for  compression  data  (solid  lines) 
and  20-pm-radius  indenter  data  (dotted  lines).  Curves  are  shown  for  (a)  Ti-6-4, 
(b)  RHA  steel,  (c)  WC-Co,  (d)  PMMA,  and  (e)  PC  (continued). 
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True  Strain  (mm/mm) 

(e) 


Figure  2.  Curves  showing  true  stress  and  true  strain  for  compression  data  (solid 

lines)  and  20-pm-radius  indenter  data  (dotted  lines).  Curves  are  shown  for 
(a)  Ti-6-4,  (b)  RHA  steel,  (c)  WC-Co,  (d)  PMMA,  and  (e)  PC  (continued). 


5.  Conclusions 


A  method  was  developed  that  allows  stress-strain  curves  to  be  detennined  from  load- 
displacement  indentation  data.  These  curves  were  shown  to  provide  a  fair  estimation  of  modulus 
and  compressive  yield  stress  information  when  compared  to  bulk  values,  with  typical  errors 
being  33%,  at  worst,  and  within  just  a  few  percent,  at  best.  Successful  measurement  was  taken 
for  materials  with  a  range  of  modulus  and  yield  properties.  The  method  developed  makes  use  of 
analyses  by  Hertz,  Tabor,  and  Francis,  by  subtracting  out  the  plastic  contribution  during  the 
indentation.  Factors  such  as  radius  of  the  sphere  and  surface  roughness  played  a  role  in  the 
adequacy  of  the  method.  For  the  materials  studied,  Ti-6-4  and  PC  had  the  greatest  agreement 
between  macroscale  and  microscale  data.  Future  work  should  include  exploring  a  more  rigorous 
framework  for  detennining  the  constraint  factor,  determining  quantitatively  the  effect  of  surface 
roughness  on  the  method’s  efficacy,  using  the  approach  for  small-scale  features,  and  studying  the 
effect  of  loading  rate  on  results. 
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